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ABSTRACT 

The  Solar  Cosmic  Ray  Ground-Level  Enhancements  (GLEs)  observed  on  11  and  15 
June  1991  were  distinctly  different  in  character.  The  small  GLE  on  11  June  was  mildly 
anisotropic  with  an  approximately  2-to-l  ratio  in  the  relativistic  proton  flux  observed  by 
"forward  viewing"  high  latitude  neutron  monitors  as  compared  with  the  flux  observed  by 
"reverse  viewing"  high  latitude  neutron  monitors.  In  contrast  the  15  June  GLE  was  almost 
isotropic  in  spite  of  the  fact  that  the  source  solar  flare  position  was  at  lHiolongitudes  that 
were  presumably  "well-connected"  to  the  earth  via  the  average  interplanetary  magnetic  field 
topology.  A  differential  power  law  in  rigidity  seems  to  fit  die  data  in  die  region  between  1 
and  6  GV  for  both  events.  For  the  11  June  GLE  maximum  our  derived  slope  is  -5.5.  For 
the  15  June  GLE  maximum  our  derived  slope  is  -6.  It  is  our  opinion  that  the  lack  of  ob¬ 
served  flux  anisotropy  during  die  15  June  GLE  is  probably  due  to  the  very  disturbed  inter¬ 
planetary  propagation  conditions  rather  than  solar  source  characteristics. 

I.  INTRODUCTION 

The  episode  of  solar  activity  that  occurred  during  June  1991  generated  a  number  of  en¬ 
ergetic  phenomena  including  intense  X-ray  and  gamma  ray  emission,  solar  neutron  emis¬ 
sion,  and  the  acceleration  of  ions  to  relativistic  energies  detectable  at  the  earth.  The  details 
of  these  solar  flares  and  their  energetic  X-ray,  gamma  ray  and  solar  neutron  emission  are 
described  elsewhere  in  this  volume  and  will  not  be  repeated  here.  The  powerful  solar  flares 
in  this  activity  episode  generated  interplanetary  shocks  that  propagated  through  the  helio¬ 
sphere.  Six  sudden  commencement  geomagnetic  storm  onsets  were  recorded  at  the  earth 
between  4  and  12  June.  The  intense  solar  activity  contributed  to  the  historic  cosmic  ray 
intensity  minimum  observed  during  this  month.  While  the  effects  of  this  activity  episode 
on  the  propagation  conditions  in  die  heliosphere  have  not  been  fully  ascertained,  all  of  these 
effects  strongly  suggest  that  propagation  conditions  were  not  quiescent.  Also  suggestive  of 
the  non-quiescent  propagation  conditions  were  the  variations  in  die  pre-event  cosmic  ray 
background  which  exceeded  the  variations  expected  from  Poisson  statistics. 

The  energy  source  for  the  solar  cosmic  ray  ground-level  enhancements  (GLEs)  studied 
in  this  paper  were  energetic  solar  flares  in  NOAA  region  6659.  The  11  June  1991  GLE  is 
time-associated  with  the  X12/3B  solar  flare  at  heliographic  coordinates  N31,  W17  having 
an  H-alpha  onset  of  0156  UT.  The  15  June  1991  GLE  is  time-associated  with  the  X12/3B 
solar  flare  at  heliographic  coordinates  N33,  W69  having  an  H-alpha  oaset  of  0810  UT. 
The  historic  cosmic  ray  intensity  low  was  recorded  on  13  June  1991,  in  the  interval  be¬ 
tween  these  two  GLEs. 
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H.  METHOD  OF  DETERMINING  HIGH  ENERGY  SOLAR  PARTICLE  SPECTRA 
FROM  THE  ANALYSIS  OF  NEUTRON  MONITOR  DATA 

An  introduction  to  the  general  concept  of  using  cosmic  ray  neutron  monitor  data  for 
the  analysis  of  high  energy  solar  proton  events  is  given  in  this  volume  by  Debrunner1  and 
that  material  will  not  be  repeated  in  this  paper.  We  have  used  our  standard  technique  for 
the  analysis  of  GLEs2  to  determine  the  spectral  characteristics  and  flux  anisotropy  of  the 
June  GLEs.  The  method  is  designed  to  reproduce  the  increase  observed  by  the  individual 
neutron  monitors  around  foe  world.  This  is  done  from  a  numerical  analysis  of  foe  solar 
particle  spectrum,  foe  flux  anisotropy,  foe  asymptotic  cone  of  acceptance  for  each  station 
and  foe  neutron  monitor  yield  function3.  For  this  analysis  we  have  used  foe  Debrunner4  et 
al.  neutron  monitor  yield  functions  to  reproduce  foe  observed  increases  at  each  neutron 
monitor.  We  have  limited  foe  form  of  foe  spectral  parameters  used  to  a  differential  power 
law  in  rigidity  which  seems  to  produce  a  satisfactory  fit  to  foe  neutron  monitor  data  in  foe 
rigidity  range  between  1  and  6  GV. 

We  model  foe  increase  utilizing  foe  functional  form, 

OO 

I  =  E  Ja(a,R)  S(R)  G(a)  AR  (1) 

whore  I  is  foe  increase  at  foe  neutron  monitor,  Rc  is  foe  cutoff  rigidity,  Jtt(a,R)  is  foe  dif¬ 
ferential  flux  in  foe  interplanetary  medium  at  pitch  angle  a  and  rigidity  R  that  is  allowed 
through  foe  asymptotic  cone  of  acceptance,  S(R)  is  foe  neutron  monitor  specific  yield  as  a 
function  of  rigidity,  and  G(cr )  is  foe  anisotropic  pitch  angle  distribution.  In  our  modeling 
approach  we  sum  foe  spectrum-yield  response  for  each  station  from  0.7  (or  foe  cutoff  rigid¬ 
ity)  at  0.1  GV  intervals  to  25  GV. 

In  describing  our  method  it  is  necessary  to  explain  foe  concept  of  asymptotic  directions 
of  approach.  Charged  particles  of  a  specified  energy  arriving  at  a  detector  from  a  specific 
direction  can  be  "mapped’  through  foe  geomagnetic  field  to  a  specific  direction  in 
space5-6’7.  The  asymptotic  direction  of  approach  defines  an  allowed  particle’s  direction  in 
space  prior  to  its  interaction  with  foe  earth’s  magnetic  field.  From  foe  "geomagnetic  op¬ 
tics"  of  high  latitude  neutron  monitors,  we  can  determine  foe  orientation  of  foe  asymptotic 
cone  of  acceptance  to  foe  interplanetary  magnetic  field  direction  and  estimate  foe  flux  ar¬ 
riving  at  each  station.  In  our  model  calculations  we  define  pitch  angle  zero  as  foe  direction 
of  foe  maximum  particle  flux  which  generally  corresponds  to  foe  direction  of  foe  inter¬ 
planetary  magnetic  field. 

DI.  THE  GLE  OF  11  JUNE  1992 

The  world-wide  network  of  neutron  monitors  on  foe  earth  recorded  a  small,  mildly 
anisotropic  GLE  on  11  June.  For  this  event  there  was  an  impulsive  onset  at  foe  forward 
viewing  stations  (those  having  asymptotic  directions  of  approach  viewing  into  foe  solar 
particle  flux  propagating  along  foe  interplanetary  magnetic  field  direction  away  from  foe 
sun)  in  foe  five-minute  interval  0235-0240  UT8.  The  onset  for  reverse  viewing  stations 
(those  having  asymptotic  directions  of  approach  viewing  into  foe  particle  flux  propagating 
along  foe  interplanetary  magnetic  field  back  toward  foe  sun)  was  after  0300  UT.  At  foe 
time  of  foe  GLE  maximum  at  about  0330  UT  stations  viewing  in  foe  probable  forward  di¬ 
rection  such  as  Apatity,  Russia  and  Mawson,  Antarctica  recorded  an  increase  of  -7  per¬ 
cent  while  stations  viewing  in  foe  probable  reverse  direction  such  as  Tixie  Bay,  Russia  and 
Inuvik,  Canada  recorded  an  increase  of  -3  percent  An  overall  conceptual  view  of  this 
small  GLE  can  be  obtained  from  Figure  1.  In  this  composite  figure  we  show  foe  asymp¬ 
totic  viewing  directions  for  selected  high  latitude  neutron  monitors  on  foe  right  The  top 
left  shows  foe  increase  observed  by  a  station  viewing  into  foe  forward  propagating  flux  and 
a  station  viewing  into  foe  reverse  particle  flux.  The  bottom  left  illustrates  foe  pitch  angle 
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distribution  required  to  generate  die  observed  particle  anisotropy  at  the  GLE  maximum8. 
This  form  is  similar  to  the  exponential  form  derived  by  Beeck  and  Wibberenz9  . 


FIG.  1.  The  11  June  1991  GLE.  Left  top :  relative  increase  observed  by  Apatity,  Russia 
(forward  viewing)  and  Tide  Bay,  Russia  (reverse  viewing).  Right:  Display  of  the  asymp¬ 
totic  viewing  directions  responsible  for  10  to  90  percent  of  the  response  of  selected  high 
latitude  neutron  monitors.  (Five  minute  data  are  displayed  until  maximum ;  15  minute  av¬ 
erages  thereafter.)  At  die  time  of  the  GLE  maximum  (0330  UT),  the  subsolar  point  was  at 
23°  N,  127°  E.  Bottom  left:  Solar  particle  flux  pitch  angle  distribution  necessary  to  produce 
the  observed  variation  in  the  high  latitude  neutron  monitors  at  the  GLE  maximum. 

Unfortunately,  there  are  no  interplanetary  magnetic  held  (IMF)  measurements  by 
earth-orbiting  spacecraft  for  11  June  1991.  However,  we  can  use  the  anisotropy  of  the  ob¬ 
served  increase  and  die  onset  time  to  approximate  die  probable  IMF  direction,  at  least  to 
the  proper  octant  The  geomagnetic  held  was  severely  disturbed,  and  a  geomagnetic  storm 
was  in  progress.  At  the  GLE  onset  the  Dst  was  -96  nT  and  increasing  toward  die  maxi¬ 
mum  of  -140  nT  which  was  observed  at  06  hours  UT.  In  our  analysis  of  this  event,  the 
IMF  direction  did  not  appear  to  be  stable  during  the  GLE,  but  in  our  opinion,  in  this  case 
this  is  not  a  serious  impediment  to  making  a  useful  spectral  determination. 

A  power  law  in  rigidity  having  a  slope  of  -5.5  yielded  a  satisfactory  fit  between  the  in¬ 
creases  of  the  various  neutron  monitors  as  a  function  of  latitude.  For  this  event  there  were 
no  significant  increases  reported  for  stations  whose  geomagnetic  cutoff  exceeded  4  GY.  A 
harder  spectrum  would  predict  increases  at  stations  having  a  cutoff  rigidity  >  4  GV.  We 
determined  the  magnitude  of  the  particle  flux  parallel  to  the  IMF  direction  and  the  flux  av¬ 
eraged  over  all  directions.  For  this  event  die  differential  power  law  in  rigidity  that  fits  the 
neutron  monitor  data  in  the  range  of  1  to  4  GV  at  the  0330  UT  GLE  maximum  is: 

J(|  =  4.93  P5-5;  and  J(lvg)  =3.49  P5-5.  (2) 

Jll  is  the  flux  in  units  of  (cm2-s-ster-GV)'1  parallel  to  the  interplanetary  magnetic  field  (i-e- 
pitch  angle  of  zero).  is  obtained  by  summing  the  anisotropic  flux  over  4r  steradians. 
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A.  Comparison  of  the  11  June  1991  GLE  Spectra  with  Spacecraft  Data. 

We  can  compare  the  spectrum  of  the  more  rigid  particles  (>  1  GV  or  >  433  MeV)  de¬ 
rived  from  the  analysis  of  neutron  monitor  data  with  the  GOES  6  and  7  five  minute  data 
We  have  integrated  the  GLE  spectrum  derived  from  the  analysis  of  the  neutron  monitor 
data  and  extended  it  to  30  MeV  as  illustrated  in  Figure  2.  The  heavy  line  in  this  figure  in¬ 
dicates  die  spectrum  that  is  derived  from  the  high  energy  flux  observed  by  the  neutron 
monitors  at  the  0330  UT  GLE  maximum  The  light  line  is  an  extension  of  this  spectrum  to 
the  spacecraft  measurement  energies.  The  +  symbol  identifies  the  spacecraft  observed 
fluxes  at  the  time  of  the  GLE  maximum 

In  addition  to  data  at  30,  SO,  60  and  100  MeV,  there  is  a  higher  energy  particle  detec¬ 
tor  on  the  GOES-6  spacecraft  from  which  integral  flux  above  355,  433  and  505  MeV  can  be 
obtained10.  We  have  plotted  these  data  in  Figure  2  for  comparison.  The  empirical  correc¬ 
tion  for  side  penetration  of  the  sensor  by  high  energy  particles  provided  by  Sauer10  has 
been  applied  to  these  data.  The  spacecraft  data  show  a  velocity  dispersive  time-of-maxima 
in  the  initial  part  of  the  event.  The  velocity  dispersive  flux  maximum  for  energies  >  300 
MeV  to  >30  MeV  occurred  after  the  0330  UT  GLE  flux  maximum,  between  0430  and 
0500  UT.  These  data  (indicated  by  die  •  symbol  in  Figure  2)  can  be  used  to  construct  a 
time-of-maximum  spectrum. 

The  GOES  spacecraft  particle  flux  data10  for  this  event  are  shown  in  Figure  3.  To  our 
biased  observations,  there  are  two  maxima  displayed  in  these  riara  a  velocity  dispersive 
flux  maximum  occurs  between  0330  UT  and  0500  UT.  There  is  a  second,  larger  and  non¬ 
velocity-dispersive  flux  maximum  that  occurs  at  about  14  hours  UT  during  an  extended  pe¬ 
riod  of  increased  magnetic  activity  when  Dst  exceeds  -100  nT.  It  is  our  opinion  that  this 
strongly  indicates  an  interplanetary  source  of  die  particles  contributing  to  the  second  maxi¬ 
mum.  In  view  of  this  we  are  reluctant  to  take  the  integrated  fluence  observed  by  earth-or¬ 
biting  satellites  for  this  event  and  extrapolate  it  back  to  the  sun  to  estimate  the  number  of 
protons  released  from  the  acceleration  site. 
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FIG.  2.  Integral  energy  spectra  for  the  GLE  of  11  June  1991.  The  heavy  dark  line  indi¬ 
cates  the  spectrum  derived  from  neutron  monitor  data.  The  light  line  is  the  extension  to 
satellite  measurement  energies.  The  +  symbol  indicates  the  spacecraft  measured  integral 
flux  at  the  time  of  die  GLE  maximum  (0330  UT).  The  •  symbol  indicates  the  spacecraft 
measured  integral  flux  during  the  0430-0500  UT  maximum. 
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FIG.  3.  The  solar  particle  flux  observed  by  the  GOES  spacecraft  for  the  11  June  1991 
event.  Note  the  non-velocity-dispersive  maximum  at  about  14  hours  UT  which  corresponds 
in  time  to  the  maximum  of  die  geomagnetic  storm. 

IV.  THE  GLE  OF  15  JUNE  1991 


The  world-wide  network  of  neutron  monitors  recorded  a  small,  approximately 
isotropic,  long-duration  GLE  on  15  June.  The  disturbed  propagation  conditions  make  it 
difficult  to  determine  precisely  die  onset  of  this  relatively  slow-rising  GLE.  The  increase 
systematically  equaled  or  exceeded  the  pre-event  background  variations  in  the  0835-0840 
UT  time  interval12.  All  high  latitude  stations  definitely  exceeded  the  pre-event  background 
variations  after  0840  UT.  We  cannot  identify  a  definite  anisotropy  in  the  onsets  of  the  for¬ 
ward  viewing  stations  as  compared  to  reverse  viewing  stations.  The  IMP-8  spacecraft 
recorded  a  velocity-dispersive  onset  in  its  measurement  range  of  8  to  400  MeV13.  The  on¬ 
set  for  the  highest  energy  measurement  (190-400  MeV  )  is  in  the  0835-0840  time  interval, 
essentially  in  time  coincidence  with  the  neutron  monitor  onsets. 

At  die  time  of  die  GLE  maximum  at  about  0930  UT  ail  high  latitude  neutron  monitors 
recorded  an  increase  of  -20  +4%.  There  was  a  very  small  flux  amplitude  anisotropy  with 
stations  viewing  in  the  forward  flux  propagation  direction  such  as  Goose  Bay,  Canada 
recording  an  —22  percent  increase  while  stations  viewing  in  the  reverse  flux  propagation 
direction  such  as  Tixie  Bay,  Russia  observed  an  increase  of  - 17  percent.  An  overall  view 
of  this  long  lasting,  approximately  isotropic  GLE  is  presented  in  Figure  4. 

There  are  direct  interplanetary  magnetic  field  (IMF)  measurements  by  earth-orbiting 
spacecraft  for  15  June  1991  until  09  hours  UT.  Then  there  is  a  4  hour  data  gap,  one  hour 
of  IMF  data  at  14  UT,  a  one  hour  data  gap,  and  then  IMF  data  are  present  for  UT  hours  16 
through  20.  Fortunately  there  did  not  appear  to  be  large  variations  in  the  hourly  averaged 
IMF  direction  and  we  assumed  that  the  IMF  direction  does  not  have  large  deviations  during 
the  data  gaps.  At  09  hours  UT  the  observed  IMF  was  at  GSE  latitude  -23°,  and  GSE  lon¬ 
gitude  145°.  Therefore  the  probable  viewing  direction  into  the  solar  proton  flux  was  at 
GSE  longitude  of  -35°. 
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FIG.  4.  Illustration  of  the  essentially  isotropic  GLE  of  IS  June  1991.  Left:  Relative  in¬ 
crease  observed  by  selected  neutron  monitors.  Right:  Display  of  the  asymptotic  viewing  di¬ 
rections  responsible  for  10  to  90  percent  of  the  response  of  selected  high  latitude  neutron 
monitors.  At  die  GLE  maximum  (0930  UT),  the  subsolar  point  was  at  23°  N,  37.5°  E. 

The  geomagnetic  field  had  been  severely  disturbed,  and  was  undergoing  a  slow  recov¬ 
ery  from  a  major  geomagnetic  storm.  At  the  GLE  onset  the  Dst  was  -41  nT  and  was  slowly 
recovering  during  the  remainder  of  the  day. 

A  power  law  in  rigidity  having  a  slope  of  -6.0  yields  a  satisfactory  fit  to  die  increases 
observed  by  die  various  neutron  monitors  as  a  function  of  latitude.  For  this  event  there 
were  measurable  increases  for  stations  at  a  quiescent  geomagnetic  cutoff  of  —6  GV.  In  our 
analysis  method,  if  we  assume  an  omnidirectional  flux,  a  differential  power  law  in  rigidity 
with  a  slope  of  -6.0  generates  the  observed  0.7%  increase  at  Rome  ^  =  6.3  GV)  and  an 
equivalent  increase  when  corrected  to  sea  level  at  die  18-NM-64  neutron  monitor  (3340 
meters  altitude)  at  Alma-Ata,  Kazakhstan  (R^  =  6.6  GV).  When  we  include  the  slight 
anisotropy  we  And  that  the  spectrum  cannot  be  harder  than  -5.5  or  we  would  predict  a 
larger  increase  than  was  observed  at  these  stations.  Since  this  was  a  long  duration  GLE  we 
can  also  determine  die  spectrum  at  other  times.  We  derive  a  differentia]  rigidity  spectrum 
which  gives  the  flux,  J,  in  units  of  (cn^-s-ster-G V)'1 .  We  find  die  high  energy  solar  cos¬ 
mic  ray  differential  rigidity  spectrum  to  be 

J  =  19.7  P-6  0  at  0930  UT,  and  J=  12.5  P*-°  at  1030  UT. 

A.  Comparison  of  the  15  June  1991  GLE  Spectra  with  Spacecraft  Data. 

In  Figure  5  we  have  integrated  the  spectrum  derived  from  die  analysis  of  the  neutron 
monitor  data  and  extended  this  spectrum  to  die  lower  energies  for  comparison  with  data 
from  the  GOES  spacecraft  The  heavy  line  in  each  panel  indicates  die  spectrum  derived 
from  die  high  energy  flux  observed  by  the  neutron  monitors;  the  left  panel  shows  the  spec¬ 
trum  at  the  0930  UT  GLE  maximum,  the  right  panel  one  hour  later.  The  light  line  is  an 
extension  of  this  spectrum  to  the  spacecraft  measurement  energies.  The  +  symbol  identi¬ 
fies  the  spacecraft  observed  fluxes  at  each  time  for  energies  of  >505  MeV,  >433  MeV 
>  355  MeV,  >  100  MeV,  >  60  MeV,  >  50  MeV  and  >  30  MeV.  The  cosmic  ray  intensity 
was  rapidly  recovering  from  its  historic  intensity  low  on  13  June  and  overwhelmed  the 
remnant  of  the  high  energy  solar  cosmic  ray  flux  after  17  June. 

Inspection  of  this  figure  shows  that  the  power  law  in  rigidity  with  a  slope  of  -6.0  does 
not  smoothly  extrapolate  to  the  spacecraft  energies  below  100  MeV.  Further  inspection  of 
this  figure  suggests  that  a  "broken  power  law*  type  of  spectrum  may  be  a  better  representa¬ 
tion  of  the  solar  particle  flux  which  evolves  extremely  slowly  for  a  "well  connected’  solar 
particle  event. 
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FIG.  5.  Spectrum  derived  for  the  GLE  of  15  June  1991  ax  0930  UT  (left)  and  1030  UT 
(right).  The  differential  rigidity  spectrum  has  been  integrated  and  converted  to  an  integral 
energy  spectrum  for  comparison  purposes.  The  heavy  dark  line  indicates  the  spectrum  de¬ 
rived  from  neutron  monitor  data.  The  light  line  is  this  spectrum  extended  to  the  satellite 
measurement  energies.  The  +  symbol  indicates  the  measured  satellite  integral  flux  at  each 
indicated  time. 
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FIG.  6.  The  solar  particle  flux  observed  by  the  GOES  spacecraft  for  the  15  June  1991  GLE 
and  associated  solar  proton  event.  The  empirical  correction  for  side  penetration  of  the  sen¬ 
sor  by  high  energy  particles  provided  by  Sauer10  has  been  applied  to  these  data. 
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The  GOES  spacecraft  particle  flux  data  for  this  event  are  shown  in  Figure  6.  In  spite 
of  the  velocity  dispersive  onset  observed  by  the  IMP-8  spacecraft13,  the  flux  maximum  has 
only  a  weak  velocity  dispersion.  The  GLE  maximum  (particles  with  energies  >  - 1  GeV) 
occurred  at  about  0930  UT,  and  the  maximum  for  energies  between  >  SO  MeV  to  >  100 
MeV  occurred  at  about  1000  UT.  It  is  our  opinion  that  the  intensity-time  profile  indicates 
that  interplanetary  diffusion  controls  both  the  spectral  evolution  and  the  intensity-time  pro¬ 
file.  Under  these  circumstances,  these  observations  at  1  AU  may  not  be  representative  of 
the  particle  source  release  profile. 
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